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We demonstrate how correlation functions for non-diagonal operators can be measured with the loop-cluster 
algorithm for quantum spin models. We introduce the U(l) quantum link model and present the construction of 
a cluster algorithm for the model. We further show how an improved estimator for Wilson loops can be realized 
with the flux-cluster algorithm. 



1. INTRODUCTION 

Cluster algorithms have been established as 
very efficient tools for the numerical simulation 
of both classical and quantum spin systems. In 
particular, for quantum spin systems, the loop- 
cluster algorithm has allowed very efficient 
studies of the 2-d Heisenberg antiferromagnet uti- 
lizing a continuous Euclidean time formulation 
j2|. An interesting development is the observa- 
tion ||] that non-diagonal correlation functions 
can be measured with the loop-cluster algorithm. 
The methodology is carried over to the numeri- 
cal study of the 17(1) quantum link model (QLM) 
which is a new, discrete formulation of Wilson's 
gauge theory Q. A flux-cluster algorithm has 
been constructed and an improved estimator for 
Wilson loops is shown to exist. 

2. THE LOOP-CLUSTER ALGORITHM 

Consider the 2-d quantum XY model, defined 
by the Hamiltonian 

H = -J ^2 i S l S x+(i + s l s l+[i\ > (!) 

x,fi— 1,2 

which couples nearest-neighbor quantum spin op- 
erators S x . The quantum partition function 

Z = Trexp(-/3H) (2) 
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is formulated as a (2+l)-d path integral with the 
system evolving in Euclidean time with extent (3. 
Time is discretized in small steps of size e and 
the Hamiltonian is decomposed in an even-odd 
pattern |2| such that the transfer matrix between 
two subsequent time slices is a product of mutu- 
ally commuting elementary transfer matrices 

T = exp( e J[S^ +A + SX +/l ])- (3) 

A complete set of states is inserted in each time- 
slice to form the path integral. For s=l/2 spins, 
the 4x4 matrix elements of T between the eigen- 
states — for example in the S 3 basis — of the two 
spins, define the Boltzmann weights for the inter- 
action between two pairs of spin states in adjacent 
time slices. The loop-cluster algorithm starts on 
a random spin state and examines the weight of 
the interacting quartette of spins. It joins proba- 
bilistically the spin to one of the other spins and 
the process continues until a loop of spin states 
is formed. The loop is the worldline of evolu- 
tion of a spin state in the (2+l)-d periodic vol- 
ume. The cluster rules are ergodic and obey de- 
tailed balance. The spin states in the cluster are 
then flipped resulting in a very efficient update 
of the system. It is shown in Q that the quan- 
tum partition function is mapped to a quantum 
random cluster model similarly to the Fortuin- 
Kasteleyn mapping for the Potts model. The 
mapping holds at the operator level and therefore 
the loop-cluster has geometrical properties inde- 
pendent of the chosen basis. Due to the discrete- 
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ness of the Hilbert space, the continuum limit 
e — > can be taken and the algorithm can be 
implemented in the continuum of the Euclidean 
time direction g]. 

The expectation value of a Green's function like 

(SlSl) = ±1r[SlSl eM-PH)} , (4) 

appears impossible to measure from the ensemble 
of configurations that contribute to Z . This is be- 
cause the insertion of the non-diagonal operator 
S 1 creates a configuration that never contributes 
to the partition function. The important observa- 
tion made in || is that the loop-cluster algorithm 
generates configurations that can be thought as 
contributing to (Q) also. Indeed, consider the loop 
which passes from the sites x and y on a fixed 
time-slice. We can imagine that we cut the loop 
on this time-slice and flip only half of the clus- 
ter. Since S 1 flips the S 3 eigenstates, this pro- 
cess precisely contributes to (Q). We therefore 
generate the loops with the algorithm and then 
perform cuts on all time-slices which dissect the 
loop in two pieces. If each half can be flipped 
independently, a +1 contribution is registered for 
the correlation function between the cutting sites. 
In this way, an improved estimator for the non- 
diagonal Greens's function ([I]) has been realized 
which does not suffer from any sign cancelations. 

3. THE U(l) QUANTUM LINK MODEL 

A new formulation for lattice gauge theories 
has been proposed in Q . Here we study the pure 
U(l) gauge theory on the lattice. We consider 
the quantization of the links, i.e. the promotion 
of the classical link field u x ^ = exp(iip X}IJ/ ) to a 
quantum link operator U x ^ acting on the link- 
based Hilbert space. The U(l) QLM is defined 
through the 4-d Hamiltonian 

H = - J E [U^U X+ ^UI + ^UU (5) 

which is invariant under the local t/(l) transfor- 
mations at the left and right ends of the link 

U' Xtll = Y\_exp(ia y Gy)U x ^Y\_ ex P(- ia zG z ) 

= eKp(ia x )U x ^exp(-ia x+ p,) . (6) 



The transformation law results in the following 
commutation relations between the local genera- 
tor of the symmetry G x and the link operators 

[Gy: U x ,ij\ — (&y.x &y,x-\-£t)Ux,l<l • (7) 

These relations are satisfied if we represent the 
link operators and the U(l) generator as 

U x ,fi = S xfl + iS xfl = S x)1 (8) 

Gx = y ] (^x.m — ^x—p,,fi) i 

where S x „ is a spin operator associated with a 
link, with the usual commutation relations 

The operator S 3 . represents the electric flux on 
the link and U Xttl acts as a flux raising operator. 
The Hilbert space per link is the space of the local 
SU(2) algebra representation and it is therefore 
finite. The quantum partition function 

Z = Trexp(-/3iT) (10) 

is formulated as a (4+l)-d path integral with the 
system evolving in a fifth dimension of extent (3. 
When j3 exceeds a critical value, the excitations 
of the model are in the Coulomb phase of the 5-d 
Abelian gauge theory. Due to the infinite cor- 
relation length in the 5-d Coulomb phase, if the 
extent of the fifth dimension (3 is finite, the the- 
ory will dimensionally reduce to the 4-d Abelian 
gauge theory. The phase transition of Wilson's 
1/(1) gauge theory can be studied from the growth 
of the correlation length near the critical fifth di- 
mension extent of the U(l) QLM. Even spin-1/2 
quantum links may suffice for the dimensional 
reduction — more on the phase diagram of the 
model for general spin in ||. 

4. THE FLUX-CLUSTER ALGORITHM 

A cluster algorithm has been constructed for 
the spin-1/2 17(1) QLM The fifth dimen- 

sion in (|l0j) is discretized in small steps of size e 
and a complete set of states is inserted at each 
time-slice. Each plane of the 4-d lattice Hamilto- 
nian is decomposed in a checker board fashion of 
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even and odd plaquettes. Only one type of pla- 
quettes is allowed to interact between two subse- 
quent time-slices so that the full transfer matrix 
is a product of commuting transfer matrices 

r = eM^W^u x+ ^ul +0 jJl u (ii) 

The 16 x 16 matrix elements of T define the Boltz- 
mann weights of the checker boarded cubes car- 
rying the interaction of 8 electric flux states e. x ,n,t 
in the path integral 



Z =U E c MS[e}) 



(12) 



The flux-cluster algorithm builds a cluster of links 
in the 5-d volume and updates the system by flip- 
ping the electric flux on the links which belong to 
the cluster. The cluster starts on a random link 
and probabilistically assigns some of the 8 links 
of the interaction cube in the cluster. The clus- 
ter rules are ergodic and obey detailed balance. 
The process continues through the 5-d volume un- 
til the cluster is completed. The cluster is a 2-d 
surface embedded in the 5-d volume and is the 
worldsheet of the evolution in the fifth dimension 
of open and closed electric flux strings. Open 
flux strings are allowed because the Gauss law 
is not enforced in ([!(]). Due to the discreteness 
of the flux, the algorithm can operate directly in 
the continuum limit e — > of the Euclidean time 
direction 

The order parameter for the gauge theory is 
the expectation value of Wilson loops Wc which is 
the product of quantum link operators U Xjfl along 
some closed curve C in the 4-d lattice at a fixed 
value of t. The link operators on the curve raise 
the electric flux and therefore the configurations 
which are generated by the algorithm will never 
contribute to the expectation value 

{Wc) = ~Tr[W c exp(-PH)]- (13) 

The flux-cluster algorithm provides the solution 
again. Each cluster which is generated and 
flipped contributes to Z, but it may also con- 
tribute to (Wc). The operator Wc flips the ori- 
entation of the electric flux along the curve C. If 



the cluster is cut along C and the flux is flipped 
on only one of the pieces, a +1 contribution to 
(Wc) is obtained. Notice that if the cluster is 
connected to itself by wrapping over the periodic 
fifth dimension, i.e. has the topology of a torus, 
it cannot be cut in two disjoint pieces along a cir- 
cle and does not contribute to (Wc)- With the 
flux-cluster algorithm therefore, a lot of informa- 
tion can be collected by examining the topology 
of the cluster and cutting the cluster along fixed 
time-slices. Each closed curve cut which allows 
independent flips of the two parts offers a +1 
contribution to the corresponding Wilson loop. 
We have therefore realized an improved estimator 
with no sign cancelations for the order parameter 
of the Abelian gauge theory. 

In conclusion, the quantum link formulation 
and the cluster algorithm have provided a com- 
plete setup for an efficient study of the Abelian 
lattice gauge theory. 
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